Determining global mean-first-passage time of random walks on Vicsek fractals 

using eigenvalues of Laplacian matrices 
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The family of Vicsek fractals is one of the most important and frequently studied regular fractal 
classes, and it is of considerable interest to understand the dynamical processes on this treelike 
fractal family. In this paper, we investigate discrete random walks on the Vicsek fractals, with 
the aim to obtain the exact solutions to the global mean-first-passage time (GMFPT), defined 
as the average of first-passage time (FPT) between two nodes over the whole family of fractals. 
Based on the known connections between FPTs, effective resistance, and the eigenvalues of graph 
Laplacian, we determine implicitly the GMFPT of the Vicsek fractals, which is corroborated by 
numerical results. The obtained closed-form solution shows that the GMFPT approximately grows 
as a power-law function with system size (number of all nodes), with the exponent lies between 1 
and 2. We then provide both the upper bound and lower bound for GMFPT of general trees, and 
show that the leading behavior of the upper bound is the square of system size and the dominating 
scaling of the lower bound varies linearly with system size. We also show that the upper bound 
can be achieved in linear chains and the lower bound can be reached in star graphs. This study 
provides a comprehensive understanding of random walks on the Vicsek fractals and general treelike 
networks. 

PACS numbers: 05.40.Fb, 61.43. Hv, 89.75.Hc, 05.60.Cd 



I. INTRODUCTION 



Fractals are an important concept characterizing the 
features of real systems, because they can model a broad 
range of objects in nature and society pj. Over the past 
few decades, fractals have attracted considerable inter- 
est from the physics community 0, Q . Among numer- 
ous fractal classes, the so-called regular fractals are an 
integral family of fractals. Examples include the Sierpin- 
ski gasket Q, the Koch snowflake the Vicsek frac- 
tals [fj|, and so on. These structures have received much 
attention , and continue to be an active object of 
research [7[. One of the main reasons for studying reg- 
ular fractals is that one can obtain explicit closed-form 
solutions on a finite structure. Another justification is 
that various problems intractable on Euclidean lattices 
become solvable on regular fractals [8|. On the other 
hand, the exact solutions on regular fractals can provide 
useful insight different from that given by the approxi- 
mate solutions for random fractals. 

A central issue, still debated, is to understand how 
the underlying geometrical and structural features in- 
fluence various dynamics defined on complex systems, 
which has been considered to be an important problem 
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in many interdisciplinary fields, e.g. network science 0- 
[Tll | . Amongst a plethora of fundamental dynamical pro- 
cesses, random walks are crucial to a lot of branches of 
sciences and engineering and have appealed much inter- 
est 0, fl^ - fl4T | . A basic quantity relevant to random walks 
is first-passage time (FPT) [15}, which is the expected 
time to hit a target node for the first time for a walker 
staring from a source node. It is a quantitative indicator 
to characterize the transport efficiency, and carries much 
information of random walks since many other quantities 
can be expressed in terms of it. Thus, a growing num- 
ber of studies have been concentrated on this interesting 
quantity (HQ. 

In view of the significance of regular fractals and 
random- walk dynamics, many authors have devoted their 
endeavors to study random walk on regular fractals I24| , 
such as the Sierpinski gasket |25l . [26j . the T~ fractal [271 - 
[2fjj , the Vicsek fractals [30l l3l|. as well as the hierarchical 
lattice fractals 32, 33]. The results of these investigations 
unveiled many unusual and exotic phenomena of random 
walks on regular fractals. But in the aspect of FPT, 
these studies only addressed the mean of FPTs between 
part of the node pairs, e.g., b etween a given node and 
all other nodes [25|, [26|, |20L f33lj , while the scaling for the 
FPT averaged over all pairs of nodes, often called global 
mean first-passage time (GMFPT), in the regular frac- 
tals is still not well understood [20j , in spite that GMFPT 
provides comprehensive information of random walks on 
fractals and other media. 

In this paper, we study analytically the discrete ran- 
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dom walks on a class of treelike fractals — Vicsek frac- 
tals, which are typical candidates for exact mathemati- 
cal fractals and have received extensive interest [33 - l38l |. 
We determine exactly the GMFPT between two nodes 
over the whole fractal family, which is verified by numer- 
ical results. The closed-form formula for the GMFPT is 
achieved iteratively by using the advantage of the specific 
construction of the Vicsek fractals. The obtained explicit 
expression indicates that for large systems the GMFPT 
increases algebraically with the size of the systems. In 
the second part of this work, we provide the rigorous up- 
per and lower bounds for GMFPT as a function of system 
size for general treelike media. We show that of all trees 
linear chains have the largest value of GMFPT and the 
star graphs have the smallest GMFPT. 



II. BRIEF INTRODUCTION TO THE VICSEK 
FRACTALS 

The so-called Vicsek fractals are constructed in an it- 
erative way d, [35J. Let V/ i9 (/ > 2, g > 1) denote 
the Vicsek fractals after g iterations (generations). The 
construction starts from (g — 1) a star-like cluster con- 
sisting of / + 1 nodes arranged in a cross-wise pattern, 
where / peripheral nodes are connected to a central node. 
This corresponds to Vf t %. For g > 2, Vf. g is obtained 
from Vf t3 -\. To obtain Vjr^, we generate / replicas 
of Vf,i and arrange them around the periphery of the 
original V/,i, then we connect the central structure by 
/ additional links to the corner copy structure. These 
replication and connection steps are repeated infinitely, 
with the needed Vicsek fractals obtained in the limit 
g — > 00, whose fractal dimension is ln(/ + 1)/ In 3. In 
Fig. [T] we show schematically the structure of V^s . Ac- 
cording to the construction algorithm, at each step the 
number of nodes in the systems increases by a factor of 
/ + 1, thus, we can easily know that the total number 
of nodes (i.e., network order or system size) of Vf l3 is 
N g = (/ + l) 9 . Since the whole family of Vicsek fractals 
has a treelike structure, the total number of links in Vf g 
is E g = N g -l = (f + iy-l. 

III. GMFPT IN THE VICSEK FRACTALS 

After introducing the Vicsek fractals V/, s , we will con- 
tinue to study numerically and analytically random walks 
performed on them, which is the primary topic of this 
present paper. The random- walk model we study is a 
simple one. Assuming the time to be discrete, at each 
time step, the walker (or particle) jumps uniformly from 
its current location to one of its neighbors. The highly de- 
sirable quantity related to random walks is the GMFPT 
starting from a source point to a given target point, av- 
eraged over all node pairs of source and target points. 

The GMFPT can be obtained numerically but exactly 
via the pseudoinverse [39l |40| of the Laplacian matrix, 



FIG. 1: Illustration of the first several iterative processes of 
a particular Vicsek fractal 1/4,3. The open circles denote the 
starting structure Va,i- 



h g , of Vf.g- The entries L 9 j of L g are defined as follows: 
the off-diagonal element L 9 j = — 1 if the pair of nodes i 
and j are linked to each other, otherwise — 0; while 

the diagonal entry — di (degree of node i). The 
pseudoinverse (denoted by h' g ) of L ff is a variant of its 
inverse matrix and is defined to be 

where e g is the jV ff -dimensional "one" vector, i.e., e g = 

(i,V--,i) T - 

The FPT between any pair of nodes in Vf >g can be 
expressed in terms of the elements, L\? , of L g . Let Tij(g) 
stand for the FPT for random walks in V/ jS , starting from 
node i to node j. Then [4l| 

N g 

Tii (9) = £ ( L » - 4 9 - L ln + L\f) Lin , (2) 
n=l 

where L 9 nn is the nth diagonal entry of L g . Thus, the 
sum, T S um(sO, f° r FPTs between all node pairs in Vf i9 
reads as 

N g 

T sum (g) = Y,Y, T ^> ( 3 ) 

3 = l 

and the GMFPT, (T) g , is 

^A^^A^^g (4) 

Using Eqs. and (U), we can compute directly the 
GMFPT (T) g of the Vicsek fractals (see Fig. ©. From 
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FIG. 2: (Color online) Global mean-first-passage time (T) g 
as a function of the iteration g on a semilogarithmic scale for 
different parameter /. The filled symbols are the numerical 
results obtained by direct calculation from Eqs. ((2} and Q, 
while the empty symbols correspond to the exact values from 
Eq. (|18p . both of which are consistent with each other. 



Fig. [3J we can see that (T) g approximately grows ex- 
ponentially in g. In other words, (T) g is a power-law 
function of network order N g obeying the scaling as 
(T) g ~ (N g ) e since N g = (f+l) 9 . It should be mentioned 
that although the expression of Eq. (Q} seems compact, 
it requires computing the inversion of a matrix of order 
N g x Ng [sec Eq. ([TJ] , which make heavy demands on time 
and computational resources for large networks. Thus, 
one can calculate directly from Eq. Q the GMFPT only 
for the first iterations. On the other hand, by using the 
method of pseudoinverse matrix it is difficult and even 
impossible to obtain the leading behavior of the exponent 
characterizing the random walks. It is thus of signifi- 
cant practical importance to seek for a computationally 
cheaper method for computing the GMFPT. Fortunately, 
the particular construction of the Vicsek fractals and the 
connection [42|, |43| between effective resistance and the 
FPTs for random walks allow us to calculate analytically 
the GMFPT and the exponent 9 to obtain rigorous solu- 
tions. 

Below we will show how to avoid the computational 
complexity of inverting a matrix. To this end, we view 
Vf ;3 as resistor networks [44| by considering each edge to 
be a unit resistor. Let Rij(g) be the effective resistance 
between two nodes i and j in the electrical networks ob- 
tained from Vf t9 . Then, according to the relation between 
FPTs and effective resistance [421 |43| , we have 



T ij (g)+T ji (g)=2E g R ij (g). 
Therefore, Eq. can be rewritten as 

T sum (g)=E g J2J^RiM- 



(5) 



(6) 



Using the previously obtained results [45|, [46( , the sum 
term on the right-hand side of Eq. ^ denoted by -R S um(<?) 



can be recast as 

Rsum (g) = J2 E % (j) = 2JV 9 Etm' ( ? ) 

i^j j = l i=2 \ 

where X [ f ] (i = 2, . . . , N g ) are all the nonzero eigenval- 
ues of Laplacian matrix, L g , of the Vicsek fractals Vf <g . 
Then, we have 



i — 9 Ai 



(8) 



Having (T) g in terms of the sum of the reciprocal of all 
nonzero Laplacian eigenvalues, the next step is to deter- 
mine this sum. 

The determination of all eigenvalues of L g can be re- 
solved by using the real-space decimation method [13, 
[48j]. Assuming that one has the eigenvalues A^ 9 ' (A^ 9 ' 

0) at generation g, then the eigenvalues A[ 9+1 ' ) of the next 
generation g+1 can be obtained through the relation [4£ 

IH3 



A 



(s+i)qG?+i) 



3)(Al 9+1) -/-l) = A^ 



(9) 



(A< 9) 0) at 



By solving Eq. (J9|, each eigenvalue A- 9 " 1 
generation g gives rise to three new and different ones at 
generation g+1, denoted by A- 9+1 \ A- 9+1 \ and A- 9+1 \ 
respectively. Moreover, the newly generated eigenvalues 
keep the degeneracy of their ancestors. Considering that 
all the nonzero eigenvalues of V/,i are A^ = 1 (i = 
2, 3, . . . , /) and Ai 1 ^ = / + 1, one can obtain all nonzero 

eigenvalues A| 9 ' of L g by iteratively solving Eq. (JSJ) g — 1 
times. 

It should be stressed that although we can provide A^ 

in a recursive way, it is difficult to write A^ in an explicit 
formula. However, in what follows we will show that the 
recursive solution to X^ allows to obtain a closed-form 
expression for the sum of the reciprocal of all nonzero 
eigenvalues of L g , denoted by A g . By definition 



A 



1 



9 ^ Ag) 

i=2 A i 



(10) 



A main goal of the following text is to explicitly deter- 
mining this sum. 

Let fi g express the set of all the N g eigenvalues of L g , 
i.e., fL 



(5) \(S) 



{Ai , A 



, A^}, where the distinctness 
of the elements has been ignored. Notice that all these 
eigenvalues are either nondegenerate or degenerate [5lj |. 

The set of the former is denoted by Q g , while the set 
of the latter is denoted by Cl g '. That is to say, il g = 



Cig includes 0, / + 1 and other eigenvalues 



generated by the "seed" aS. x = / 

1 and other eigenvalues derived from 1. Furthermore, 



1; and il g 2 ' includes 
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the degeneracy of the degenerate eigenvalues rests with 
the generation at which they appeared at the first time. 
At a given generation j, the degeneracy of eigenvalues 1 
is Aj = (f — 2)(/ + l)^ 1 + 1, a degeneracy that their 
descendants keep. In what follows, for convenience we 
use Qg to represent the nondegenerate eigenvalues of 
Jj g other than 0. 

We now return to derive A g , which can be evidently 
recast as 



A 9 = 



E .(g) + E .(g) 



A (9) efi< 2) ^ 



We denote the two sums on the right-hand side of 
Eq. (fTTj) by Ag 1 -*, and A g 2 \ respectively. Below we will 
calculate the two quantities Ag and Ag 2 ' 1 . 

We first calculate Ag 1 -*. At the initial generation 1, 
there is only one nondegenerate eigenvalue / + 1 , which 
produces three different nondegenerate eigenvalues at 
generation 2. We call these three eigenvalues the first- 
generation descendants of / + 1, which give rise to 3 2 
second-generation descendants of / + 1 at the third gener- 
ation. Thus, at ith generation, 3 t ~ 1 (i — l)th-generation 
descendants of / + 1 are produced. Since all eigenval- 
ues (degenerate or nondegenerate) which appeared at 
one generation will still appear in all subsequent gen- 
erations [49l - [5lj |. we have C ^g 1 "*- Hence, as noted 

above, Jig 1 ' 1 consists of f + 1 and all its offspring produced 
after generation 1. 

Let rf 3 be the sum of the reciprocal of all the (i — 
l)th-generation descendants of / + 1. Then, Ag 1 -* can be 
rewritten in terms of as 



A (D 
a 



9-1 

E^ 

i=0 



(12) 



where = A^ = l/(/ + 1). 

Note that for each nonzero eigenvalue (degenerate or 

nondegenerate) € fig, Eq. © can be rewritten in an 
alternative way as 

( A 0+i))3 _ (/ + 4)(A 4 (9+1) ) 2 + 3(/ + 1)A< 9+1) - A 



O) 



0. 
(13) 



According to the Vieta's formulas, the three roots (i.e., 
^ 9 i +1) ' A ^2 +1) ' and \ (S 3 +1) ) 01 E q- <QD satisf y the fol- 



lowing two relations: Xf^ \ 2 
\0- 

\,2 



O+i) . \0+i) . A te 3 +1) 



A„ (9) and 



A 0+D . A 0+D + A 0+D . A 0+D + A 0+D . A 0+D = 3 (/ + 1). 



Thus, 1/A^ +1) + 1/A| 9 2 +1) + 1/A& +1) = 3(/ + 1)/A| S) . 
Based on the results obtained above, we have 

3(/ + 1) =3(/ + l)rW, (14) 



r(i) = 
g 



E 



>0) 



which together with the initial condition = l/(/ + 1) 
leads to Tf' — 3 9 (f + l) 9 ~ l . Inserting this result into 



1 3 9 (/ + 1)9-1 
7 + 1 3/ + 2 



(15) 



Eq. (QUI, we get 

a^EW+i)*- 1 ] 

After obtaining Ag 1 -*, all that is left to find an expres- 

(2) 

sion for A g is to evaluate A y g ' . For each eigenvalue 1, 
applying an approach similar to that used above, we 



can compute the sum of the reciprocal of its (i — l)th- 



(11) generation descendants, which we represent by T,[ 2 \ Af- 

(2) 

ter some simple algebra, we obtain = 3*(/ + 1)' 

(0 < i < g - 1), where T Q 2) = 1 express the recipro- 
cal of the "seed" eigenvalue 1 itself. It has been shown 
that (49l - l5l1 | in the Vicsek fractals Vf, g , the degeneracy 
of eigenvalues 1 is A g = (/ - 2)(/ + l)^ 1 + 1, and the 
degeneracy of each of its ith-generation (0 < i < g — 1) 
offspring is Ag_ 4 = (/ - 2)(/ + 1)9-^ + 1. Then, the 



quantity Ag is evaluated as follows: 



9-1 



,(2) 



Af = E Km • T 



= (/-2)(/ + l)g- 1 (3g-l) 3 g(/+l)g-l 

2 3/ + 2 " P> 

Plugging Eqs. (US]) and JTB]) into Eq. ((TTJ» yields 



A 3 = 



(/-2)(/ + l)f- 1 (3 9 -l) , / + 23»(/ + l)»-l 



f + 1 3/ + 2 



Using the relation (T) g — 2A 9 , we have 



(17) 



(T) g = (f-2)(f+ir-\v-i)+ 



2(/ + 2) 3»(/+l)»-l 



/ + 1 



3/ + 2 



We have confirmed this closed-form expression for (T) s 
against direct computation from Eqs. ([2]) and ((H). For 
all range of g and different /, they completely agree with 
each other, which shows that the analytical formula pro- 
vided by Eq. (fl"8|) is right. Figure [2] shows the compar- 
ison between the numerical and predicted results, with 
the latter plotted by the full expression for the sum in 
Eq. (HID. 

We can also support the validity of Eq. (TT8l) by us- 
ing another method. In fact, the correctness of Eq. (fT8|) 
depends on all the nonzero Laplacian eigenvalues, the 
exactness for derivation of which can be established ac- 
cording to the relation between the Laplacian eigenval- 
ues and the number of spanning trees of a graph. It 
has been established that the number of spanning tress 
on a connected graph G with order N, N at (G), is re- 
lated to all its nonzero Laplacian eigenvalues Xi (assum- 
ing Ai = and Aj 7^ for i = 2, • • • , N), obeying the 

relation N st (G) — JIi=2 El- Since the Vicsek frac- 
tals Vf g are trees, the product of all nonzero Laplacian 
eigenvalues for Vf. g , denoted by Q g , should equal N g , 
which can be corroborated by the following argument. 
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By definition, B g = • 9^ ; , where 9 



»(2) 



1,2) is 



the product of Laplacian eigenvalues in fi a . Applying 
the Vieta's formulae, we can easily obtained the product 
of the ith-order (0 < i < g — 1) offspring of the "seed" 
eigenvalue / + 1 is / + 1, which is independent of i. Then 



e 



Similarly, we have 91, 



1. Hence, 



Q g = (/ + l) 9 = A g , which proves the correctness of the 
computation on the Laplacian eigenvalues for Vf g . 

We proceed to show how to represent GMFPT, (T) g , 
as a function of the network order N g , with the aim to 
obtain the relation between these two quantities. Recall- 
ing N g = (f + 1)9, we have 3 9 = (A s ) ln3 / ln(/+1) that 
enables one to write (T) g in the following form: 



CO, 



y^N g [(N g ) 

2(/ + 2) 
(/ + l)(3/ + 2) 



ln3/ln(/+l) _ 



l(N 9 



a+ln3/ln(/+l) _ 



1119) 



Equation (fT9]) unveils the explicit dependence relation 
of GMFPT on network order N g and parameter /. For 
large systems, i.e., N g —¥ oo, we have following expression 
for the dominating term of {T) g : 



CO* 



/(3/ - 2) 

-(N U+lnS/lntZ+l) 

(/+l)(3/ + 2) 1 9) 



/(3/ - 2) 
(/+l)(3/ + 2) 

/(3/-2) 
(/ + l)(3/ + 2) 



(N a ) 9 



(20) 



where d = 21n(/ + l)/ln(3/ + 3) is the spectral di- 
mension of the Viskek fractals [5l|. Thus, in the large 
limit of 5, the GMFPT grows approximately as a power- 
law function of network order N g with the exponent 
8 = 1 + ln 3/ ln(/ + 1) being a decreasing function of /. It 
is easy to see that the exponent 8 is larger than 1 but not 
greater than 2. Particularly, when / = 2, 8 reduces to 2, 
which is the highest one reported thus far. In fact, 2 is 
largest exponent for GMFPT of random walks defined on 
treelike media, the rigorous proof of which will be given 
in the next section. In addition, it should be mentioned 
that the obtained superlinear dependence of GMFPT on 
the network order is in contrast with the other scalings 
previously observed for other media, e.g., linear scaling 
for the Apollonian networks 53j and the pseudofractal 
scale- free web |54[, a logarithmic correction to the lin- 
ear dependence for small- world trees [U [23[ . Figure [3] 
shows how the GMFPT scales with the network order for 
various parameter /. 



IV. BOUNDS FOR GMFPT IN TREES 

In Sect. Ill, we have shown that the GMFPT in the 
Vicsek fractals scales as a power-law function of network 
order. Previous studies exhibited that GMFPT in other 
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FIG. 3: (Color online) Global mean first-passage time (T) g 
versus the network order N g on a log-log scale. The filled 
symbols described the analytic results shown in Eq. (|18p . The 
solid lines represent the corresponding leading scaling given 
by Eq. ipDjl. 



trees may depend on network order N following different 
scalings. For example, in the deterministic uniform re- 
cursive trees, their GMFPT varies with network order N 
as N In N [55]; in the T- fractal, the GMFPT grows as 
jyi+ln2/ln3 jggj These show that in different trees, the 
GMFPT obeys different dependence relation on network 
order. Then, some natural questions arise: what are the 
upper and lower bounds for GMFPT in general trees? In 
which trees are these bounds reached? 

As a matter of fact, the above questions are equivalent 
to find the upper and lower bounds for the total effective 
resistance, i? SU m, as defined similarly by Eq. (0. One 
can prove with ease using various methods [57H6 1| that 
for trees with order N, the upper and lower bounds for 
^sum are 



Jj>Upp _ 



N(N- l)(N+ 1) 



and 



i? s L ™ = 2(N l) 2 



(21) 



(22) 



respectively. 

The upper bound can be only reached for the tree that 
is exactly a linear chain (a path), which has two nodes 
with degree 1 at both ends of the chain and N — 2 nodes 
with degree 2 in the middle [60(. Actually, this linear 
chain is one of the particular Viscek fractals correspond- 
ing to / = 2. The result provided by Eq. (|21|) is com- 
patible with that of the Vicsek fractals corresponding to 
/ = 2. As for the lower bound, it can be only achieved 
when the tree is a star graph [H, [H, HlJ] , consisting of one 
central node and N—l leaf nodes. All these leaf nodes are 
linked to the central node, and there is no edge between 
the leaf nodes. 
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From Eqs. ([21]) and (I22j) . we can easily obtain that the 
upper and lower bounds for GMFPT are 



(T) 



Upp_ (N-l)(N + l) 



and 



2(JV- l) 2 
N 



(23) 



(24) 



respectively. Thus, the GMFPT (T) for general trees 
satisfies the relation (T) u pp < (T) < (T) Low . For large 
trees (i.e., N — > oo), the leading scalings for (T) Upp 
and (T) Low change separately with network order N as 
(T) u pp - TV and (T) Low - TV 2 , implying that the scaling 
for the GMFPT in any tree must lie between linear scal- 
ing and square of network order. It is very obvious that 
the upper bound is much larger than the lower bound, 
the reasons for which lie with the underlying structures 
of the corresponding graphs: the linear chain is homoge- 
neous, while the star graph is heterogeneous. 

In the star graphs, the central node has a very large 
degree, and thus plays a crucial role in keeping the whole 
graph together. When the random-walk process is per- 
formed in the star graphs, the walker has a tendency to 
migrate toward the central node, through which it jumps 
to the target nodes. Therefore, the efficiency of random 
walks is very high in the star graphs, the linear scaling of 
the GMFPT with N is the best we can see jej. Notice 
that the same scaling has been previously observed for 
complete graphs 54|. In fact, the star graphs can be ob- 
tained from the complete graph with the same order by 
whittling down complete graphs, i.e., by the judiciously 
removing edges from complete graphs leaving only one 
node with N — 1 connections, in order that the walker in 
the star graphs can find the destination nodes as easily 
as in the complete graphs. 

On the contrary, in the linear chains all nodes are ho- 
mogenous. When the walker starting from the source 
point to find the target node far away from the staring 
point, it must traverse all nodes between the starting 
point and the destination node. This makes the traverse 
time much longer than in the star graphs. 

Finally, we should stress that although the star graphs 
are extreme of heterogenous media, they are very in- 
structive to understand the dynamics of random walks 



on other heterogeneous graphs, especially scale-free net- 
works [62| . which are ubiquitous in real natural and 
social systems loll] . Previous studies have shown 
that random walks in scale-free networks are very effi- 
cient [33j, [5J, I65l469| , the roots of which is actually can be 
heuristically explained as above. The large-degree nodes 
in scale-free networks play a similar role as that of the 
central node in the star graphs, making the GMFPT very 
small. 

V. CONCLUSIONS 

We have studied the discrete random walks on the fam- 
ily of Vicsek fractals, which includes the linear chains 
as a particular case. Using the connection between the 
FPTs and the Laplacian eigenvalues for general graphs, 
we have computed the GMFPT averaged over all pairs 
of nodes in the fractals and obtained explicit solutions 
to the GMFPT. The obtained closed-form formula shows 
that in the limit of infinite network order N, the GMFPT 
(T) grows approximately as a power-law function of TV: 
(T) ~ jv 1+ln3 / ln ( /+1 ). We have also provided rigorous 
bounds on the network order dependence of the GMFPT 
in general treelike networks. We showed that the upper 
and lower bounds can be achieved in linear chains and 
star graphs, respectively. Our study sheds useful insights 
into the random-walk process occurring on treelike me- 
dia. 
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